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1  Introduction 

This  is  the  final  technical  report  on  ARO  contract  DAAH04-95- 1-0100,  which  began  March 
1,  1995  and  ended  June  30,  1998.  The  scientific  focus  of  this  project  lies  at  the  frontier 
where  mathematics  meets  materials  science.  We  are  concerned  with  the  effective  moduli 
of  composites,  the  formation  of  microstructure  in  coherent  phase  transitions,  and  random 
solutions  of  nonlinear  evolution  equations.  These  diverse  problems  are  linked  by  a  common 
theme:  the  relationship  between  microstructure  and  macroscopic  behavior.  The  relevant 
mathematical  tools  include  homogenization,  the  calculus  of  variations,  and  stochastic  pro¬ 
cesses. 

The  scope  of  this  effort  is  rather  broad.  One  goal  is  to  explore  the  mechanics  of  smart 
materials.  Projects  in  this  area  include  analysis  of  relaxation  effects  observed  in  certain 
shape-memory  alloys  (Bhattacharya  and  Swart,  with  James  -  Section  2.2);  simulation  of  the 
microstructural  pathway  of  a  stress-induced  phase  transformation  (Killough,  working  with 
Kohn  -  Section  2.3);  and  a  model  for  how  symmetry  and  texture  determine  the  recoverable 
strains  of  shape-memory  polycrystals  (Bhattacharya  and  Kohn  -  Section  2.4).  A  second 
goal  is  to  explore  micromagnetics,  plasticity,  and  other  types  of  material  behavior  involving 
hysteresis  and  pattern  formation.  Projects  in  this  area  include  an  analysis  of  polycrystal 
plasticity,  when  the  basic  crystal  has  few  slip  systems  (Kohn  and  Little  -  Section  2.5);  and 
an  explanation  for  the  branching  of  magnetic  domains  in  a  uniaxial  ferromagnet  (Choksi 
and  Kohn,  with  Otto  -  Section  2.7).  A  third  goal  is  the  improved  analysis  and  design  of 
technologically  significant  composite  materials.  Projects  in  this  area  include  analyses  of 
the  dielectric  properties  of  various  heterogenous  polymer  and  liquid  crystal  systems  (Levy, 
with  Stroud  and  Palffy-Muhoray  -  Section  2.6);  and  a  model  for  the  acoustic  response  of 
a  composite  when  the  wavelength  of  the  incident  radiation  is  not  well-separated  from  the 
length  scale  of  the  microstructure  (Avellaneda,  with  Berlyand  and  Clouet  -  Section  2.1). 
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The  mentoring  of  postdoctoral  scientists  is  a  major  part  of  our  activity.  This  grant 
provided  partial  support  for  postdoctoral  researchers  Pieter  Swart  (spring  1995),  Tom  Little 
(1995-6),  Ohad  Levy  (1995-6  and  1996-7),  Rustum  Choksi  (1995-6  and  1996-7),  and  Tim 
Schulze  (spring  1998).  About  80%  of  our  budget  was  spent  on  postdoc  salary,  fringe  benefits, 
and  associated  overhead.  Our  postdocs  have  done  very  well  in  the  job  market:  Swart  has 
been  a  Staff  Scientist  at  Los  Alamos  National  Laboratory  since  1995;  Tom  Little  took  a 
postdoctoral  position  at  the  University  of  Minnesota  (1996-7)  followed  by  an  NSF  Industrial 
Postdoctoral  Fellowship  at  Columbia  University  and  Morgan  Stanley  (1997-8);  Ohad  Levy 
took  a  postdoctoral  position  at  the  Levich  Institute  of  CCNY-CUNY  (1997-8);  Rustum 
Choksi  took  a  tenure-track  Assistant  Professorship  at  Simon  Fraser  University  (1997-8). 
Tim  Schulze  is  still  at  NYU,  supported  by  a  grant  from  NSF  for  one  more  year. 

We  are  also  active  in  the  training  of  students.  We  supervised  the  thesis  research  of 
Reade  Ryan  (PhD,  July  1996),  Weimin  Jin  (PhD,  August  1997)  and  Matt  Killough  (PhD, 
August  1998).  This  contract  did  not  provide  academic  year  support  for  these  students,  but 
it  did  provide  summer  support  for  Ryan  and  Jin,  and  Ryan  received  academic-year  support 
from  an  ARO-sponsored  AASERT  contract. 

For  an  interdisciplinary  effort  such  as  ours  to  succeed,  contact  with  scientists  in  the 
relevant  application  areas  is  of  paramount  importance.  Therefore  we  often  attend  and 
sometimes  organize  conferences  at  the  interface  between  mathematics  and  mechanics.  For 
example,  in  1997  Kohn  co-organized,  with  Steve  Cowin  of  City  College,  a  Workshop  on 
Tissue  Adaptation  and  Structural  Optimization  (NYU,  1/97);  he  chaired  the  organizing 
committee  for  the  Second  SIAM  Conference  on  Mathematical  Aspects  of  Materials  Science 
(Philadelphia,  5/97);  he  was  an  invited  speaker  at  the  78th  Statistical  Mechanics  Meeting 
(Rutgers,  5/97);  he  served  on  the  Research  Advisory  Board  of  the  ARO  MURI  on  Control 
and  Design  of  Smart  Structures  (Harvard,  8/97);  and  he  attended  a  meeting  on  Multiscale 
Materials  Prediction  -  Fundamentals  and  Industrial  Applications  (MIT,  9/97). 

2  Major  Research  Accomplishments 

2.1  Piezoelectric  composites 

•  M.  Avellaneda  and  P.  Swart,  “Caculating  the  performance  of  1-3  piezoelectric  com¬ 
posites  for  hydrophone  applications:  an  effective  medium  approach,”  J  Acoust.  Soc. 
Amer.  103  (1998)  1449-1467. 

•  M.  Avellaneda,  L.  Berlyand,  and  J.-F.  Clouet,  “Frequency  dependent  acoustics  of 
composites  with  interfaces,”  preprint. 

The  work  of  Avellaneda,  Berlyand,  and  Clouet  addresses  the  scattering  and  transmission 
of  an  acoustic  wave  by  a  composite  medium  consisting  of  a  periodic  array  of  rods  embedded 
in  a  homogeneous  matrix.  The  motivation  for  this  problem  lies  in  the  analysis  and  design 
of  piezoelectric  composites  for  use  as  acoustic  transducers.  Let  A  be  the  wavelength  of  the 
incident  radiation,  and  d  the  length  scale  of  the  composite.  The  limit  d  •C  A  is  standard: 
effective  medium  theory  can  be  used,  i.e.  the  composite  can  be  treated  as  a  homogeneous 
effective  medium.  In  applications,  however,  it  is  often  important  to  consider  corrections  due 
to  nonzero  d/X.  This  is  the  “transitional”  regime,  when  frequency  dependent  effects  are  no 
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longer  negligible  (as  in  the  quasistatic  case)  but  d/X  is  small,  so  that  the  incident  frequency 
is  below  the  “cross-talk”  regime  in  which  the  rods  interact  strongly  (d  ~  A).  The  work  of 
Avellaneda  et  al.  determines  the  acoustic  properties  of  the  composite  in  the  transitional 
regime,  by  means  of  an  expansion  in  powers  of  e  =  d/X.  Each  term  in  the  expansion  is 
determined  by  solving  a  sort  of  “cell  problem.”  The  first  few  terms  determine  the  essential 
frequency-dependence  of  the  composite’s  impedance  and  reflection  coefficients. 

The  article  by  Avellaneda  and  Swart  reports  work  done  under  our  previous  ARO  grant; 
its  content  was  included  in  the  final  report  on  that  grant. 

2.2  Relaxation  effects  in  shape-memory  alloys 

•  K.  Bhattacharya,  R.D.  James,  and  RJ.  Swart,  “Relaxation  in  shape-memory  alloys. 
Part  I:  Mechanical  model,”  Acta  Materialia  45  (1997)  4547-4560. 

•  K.  Bhattacharya,  R.D.  James,  and  P.J.  Swart,  “Relaxation  in  shape-memory  alloys. 
Part  II:  Thermomechanical  model  and  proposed  experiments,”  Acta  Materialia  45 
(1997)  4561-4568. 

This  work  represents  the  first  serious  attempt  to  model  the  relaxation  effects  observed 
in  certain  shape  memory  materials  in  the  martensitic  phase.  In  these  alloys,  a  twinned 
specimen  can  easily  be  detwinned  by  loading.  If  the  loads  are  released  soon  after  detwinning, 
the  specimen  spontaneously  returns  to  the  twinned  configuration.  If  the  loads  are  held  for 
some  time,  however,  the  specimen  remains  detwinned  even  after  loads  are  released.  This 
behavior  gives  rise  to  a  wide  variety  of  effects,  including  changes  in  the  hysteresis  loops 
with  cycling.  There  is  some  controversy  as  to  its  origin,  which  may  be  due  to  transient 
shuffle- twinning  or  possibly  to  short-range  ordering.  The  work  of  Bhattacharya,  James, 
and  Swart  develops  a  simple  model  based  on  the  presence  of  a  slowly-relaxing  bulk-energy 
effect.  Numerical  simulations  show  good  qualitative  agreement  with  the  experimental  data. 
Experiments  are  proposed  which  would  further  test  the  model  and  shed  light  on  the  physical 
mechanisms  of  relaxation. 

2.3  Stress-induced  twinning  and  pattern  formation 

•  R.D.  James,  R.V.  Kohn,  and  T.W.  Shield,  “Modeling  of  branched  needle  microstruc¬ 
tures  at  the  edge  of  a  martensite  laminate,”  in  Proceedings  ICOMAT-95,  Journal  de 
Physique  IV  -  Collogue  C8  (1995)  253-259. 

•  M.  Killough,  “A  Diffuse  Interface  Approach  to  the  Development  of  Microstructure  in 
Martensite,”  Ph.D.  Thesis,  New  York  University,  in  preparation. 

Hysteresis  in  shape-memory  materials  is  associated  with  locally-stable  microstructures 
and  transitions  between  them.  Recent  experiments  by  Chu  and  James  observed  the  stress- 
induced  phase  transformation  from  one  variant  of  CuAlNi  martensite  to  another.  They  ob¬ 
served  striking  and  remarkably  systematic  microstructural  patterns  which  form  the  pathway 
for  this  transformation. 
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The  work  of  James,  Kohn,  and  Shield  explains  how  general  features  of  the  observed 
patterns  -  e.g.  the  structure  of  the  approximate  interfaces  -  is  consistent  with  a  model 
based  on  energy  minimization. 

The  Ph.D.  research  of  Killough,  advised  by  Kohn,  takes  a  different,  more  dynamic  and 
computational  approach  to  the  same  problem.  Killough’s  simulations  are  based  on  a  Landau 
theory  (with  strain-gradient  terms)  and  viscoelastodynamics  (without  inertial  effects).  They 
capture  many  qualitative  features  of  the  observed  patterns,  including  the  nucleation  and 
splitting  of  needles.  These  simulations  have  diffuse  phase  interfaces;  Killough’s  approach 
can  be  viewed  as  a  sort  of  phase-field  model,  with  strain  as  the  order  parameter.  It  is 
natural  to  inquire  about  the  associated  shaxp-interface  dynamics.  The  answer  is  surprisingly 
classical:  the  normal  velocity  of  a  twin  boundary  is  a  specific  function  of  the  orientation 
and  “driving  force.”  Recent  work  by  Abeyaratne,  Chu,  and  James  proposed  that  the 
experimentally  observed  patterns  represented  local  minima  of  (elastic  -I-  surface)  energy. 
Killough’s  simulations  call  this  into  question,  since  his  patterns  are  similar  to  those  observed 
experimentally  but  they  are  transient  rather  than  stationary. 

2.4  Shape-memory  polycrystals 

•  K.  Bhattacharya  and  R.V.  Kohn,  “Recoverable  strains  in  shape-memory  polycrystals,” 
in  Proceedings  ICOMAT-95,  Journal  de  Physique  IV  -  Collogue  C8  (1995)  261-266. 

•  K.  Bhattacharya  and  R.V.  Kohn,  “Symmetry,  texture,  and  the  recoverable  strain  of 
shape-memory  polycrystals,”  Acta  Materialia  44  (1996)  529-542 

•  K.  Bhattacharya  and  R.V.  Kohn,  “Elastic  energy  minimization  and  the  recoverable 
strains  of  polycrystalline  shape-memory  materials,”  Arch.  Rational  Mech.  Anal.  139 
(1997)  99-180. 

•  K.  Bhattacharya,  R.V.  Kohn,  and  S.  Kozlov,  “Some  examples  of  nonlinear  homoge¬ 
nization  involving  nearly  degenerate  energies,”  to  appear  in  Proc.  Roy.  Soc.  London 
A 

•  R.V.  Kohn  and  B.  Niethammer,  “Geometrically  nonlinear  shape  memory  polycrystals 
made  from  a  two- variant  material,”  in  preparation. 

•  K.  Bhattacharya,  R.V.  Kohn,  and  Y.-C.  Shu,  “The  Taylor  estimate  of  recoverable 
strains  in  shape-memoiy  polycrystals,”  in  proceedings  of  lUTAM  Symposium  on 
Transformation  Problems  in  Composites  and  Active  Materials.,  Y.A.  Bahei-El-Din  and 
G.J.  Dvorak,  eds.,  Kluwer,  in  press. 

The  work  of  Bhattacharya  and  Kohn  is  the  culmination  of  a  project  begun  in  1991, 
when  Bhattacharya  came  to  NYU  as  an  ARO-supported  postdoc.  This  work  addresses  the 
effective  behavior  of  polycrystals  made  from  a  shape-memory  material.  The  overall  goal 
is  to  explain  how  the  symmetry  of  the  underlying  phase  transformation  and  the  texture 
of  the  polycrystal  combine  to  determine  the  recoverable  strain  of  the  polycrystal.  There 
is  an  elementary  “Taylor  estimate,”  analogous  to  the  one  introduced  by  G.I.  Taylor  in 
the  theory  of  polycrystal  plasticity.  The  paper  in  Acta  Materialia  explores  this  estimate 
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and  Its  consequences;  the  one  in  Archive  for  Rational  Mechanics  and  Analysis  explores 
the  accuracy  of  the  Taylor  estimate  by  using  techniques  from  nonlinear  homogenization 
One  practical  conclusion  is  that  even  in  NiTi,  the  most  common  shape-memory  material’ 
favorable  texture  is  crucial  for  good  shape-memory  behavior.  This  is  docume^ed  in  he 
paper  by  Bhattacharya.  Kohn.  and  Shu  (and  to  other  recent  work  by  Bhattacha^^: 

The  paper  ^th  S.  Kozlov  addresses  a  special  class  of  nonlinear  homogenization  problems 
n  the  plane.  The  rnicrostructure  is  a  sort  of  checkerboard  polycrystal  and  the  ener^  of 
the  basic  crystal  is  degenerate  in  one  direction.  The  main  results  are  a  pair  of  matSine 
upper  and  lower  bounds  for  the  homogenized  energy.  This  example  shows  that  certain 
estimates  proved  by  Bhattacharya  and  Kohn  are  nearly  sharp.  It  can  also  be  understood  as 

from  >,r  conductivity,  involving  a  checkerboard  polycrystal  made 

from  a  highly  anisotropic  nonlinear  conductor.  As  current  traverses  the  polycrystal  it 

LtspotJ.  logarithmic  spirals  centered  at  certain  corners  Ihlh  form 

The  work  of  Bhattacharya  and  Kohn  uses  geometrically  linear  elasticity.  Shape-memorv 
materials  sometimes  recover  as  much  as  10%  strain,  so  the  effects  of  geometric  nonlinearitv 

and  Niethammer  provides  the  first 
extension  of  (some  of)  this  work  to  the  geometrically  nonlinear  setting.  Mattacharya  Ld 

Kohn  considered  the  recoverable  strain  of  a  2D  polycrystal,  when  the  basic  cryXl  has 
JUS  two  variants  of  martensite;  their  conclusion  was  that  for  polycrystals  with  Efficient 

“  ‘he  analogous  uonltoear  sotting 

2.5  Polycrystal  plasticity 

*  Me  ctoZ”  P'-ohlems  otpolyerystal  plasticity  with  doUciont 

oasic  crystals,  oIAM  J.  Applied  Math.^  m  press. 

This  work  addresses  the  yield  stress  of  a  (rigid,  perfectly  plastic)  polycrystal  in  terms 

first  ext  the  yield  properties  of  the  basic  crystal.  That  is  a  classical  problem  -  the 

first  extensive  discussion  was  given  by  G.I.  Taylor  in  the  1930’s  -  however  it  is  still  not  well 
understood.  Kohn  and  Little  concentrate  on  the  case  when  the  basic  crystal  has  a  highly 
eccentric  yield  set.  The  elementary  Sachs  (constant-stress)  and  Taylor  (constant-strS 
bounds  are  very  far  apart  in  this  case,  so  it  is  natural  to  ask  whether  they  can  be  improved 
and  which  one  describes  the  “typical”  behavior.  Kohn  and  Little  show,  L  antiplZ  sW 
the^SpS^  stress  versions  of  this  problem,  that  the  Bishop-Hill-Taylor  can  be  improved  while 

compaJlssTndThT?  homogenization,  compensated 

compactness  and  the  linear  comparison  method”  for  bounding  the  behavior  of  nonlinear 

thM^thirof  behavior  is  closer  to  that  of  the  Sachs  bound 

tnan  that  of  tho  Bishop-Hill-Taylor  bound. 

2.6  Linear  and  nonlinear  composites 

.  W.  Macevoy  Jr.  and  M.  Avellaueda,  "Electoosmotic  coupUug:  incorporattog  larger 
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surface  effects  with  a  new  length  scale,”  J.  Colloid  and  Interface  Science  188  (1997) 
139-149. 

•  D.J.  Bergman  and  O.  Levy,  “Nonlinear  behavior  and  optical  bistability  in  compos¬ 
ite  media,”  in  Nonlinear  Optical  Materials,  IMA  Volumes  in  Mathematics  and  its 
Applications  101,  Springer- Verlag,  1998. 

•  0.  Levy,  “Weakly  nonlinear  conductivity  and  flicker  noise  near  percolation”  in  Math¬ 
ematics  of  Multiscale  Materials,  IMA  Volumes  in  Mathematics  and  its  Applications 
99,  Springer- Verlag,  1998. 

•  R.V.  Kohn  and  0.  Levy,  “Duality  relations  for  non-Ohmic  composites,  with  applica¬ 
tions  to  behavior  near  percolation,”  J.  Stat.  Phys.  90  (1998)  159-189. 

•  0.  Levy  and  D.  Stroud,  “Macroscopic  disorder  and  the  metal-insulator  transition  in 
conducting  polymers,”  J.  Phys.:  Condensed  Matter  9  (1997)  L599— L605. 

•  O.  Levy  and  D.  Stroud,  “Maxwell  Garnett  theory  for  mixtures  of  anisotropic  inclu¬ 
sions:  applications  to  conducting  polymers,”  Phys.  Rev.  B  56  (1997)  8035-8046. 

•  O.  Levy  and  P.  Palffy-Muhoray,  “Dielectric  properties  of  polymer  dispersed  liquid 
crystals,”  in  preparation. 

•  O.  Levy  and  P.  Palffy-Muhoray,  “Dielectric  response  of  polymer  dispersed  liquid  crys¬ 
tals  near  the  nematic-isotropic  transition,”  in  preparation. 

The  work  of  Macevoy  and  Avellaneda  addresses  electroosmosis,  the  phenomenon  whereby 
an  applied  current  produces  flow  in  a  porous  medium.  Microscopically  one  is  concerned  with 
a  porous  medium  such  as  sandstone,  whose  pore  space  is  filled  with  a  fluid  containing  ions 
such  as  salt  water.  The  governing  equations  are  Stokes  flow  in  the  fluid.  Coulomb’s  equation 
for  the  electric  potential,  and  a  difi'usion  equation  for  the  ions.  Macroscopically  one  expects 
a  law  of  the  form 


Je  =  LiiV  (f>  +  L12S/P 

J{  =  iyi2V(/i -H  L22VP 

relating  the  electric  current  Je,  the  fluid  current  Jf,  the  electric  field  V^,  and  the  pressure 
gradient  Vp.  If  I/12  =  0  this  reduces  to  Ohm’s  law  and  Darcy’s  law.  The  coefficient 

Li2  is  responsible  for  electroosmosis;  it  is  nonzero  when  the  microscopic  model  includes  a 

smrface  charge.  Crucial  physical  parameters  are  the  ion  number  density  and  the  pore  size. 
After  nondimensionalization  these  determine  the  reduced  surface  charge  ps  and  the  Debye 
parameter  6.  The  work  of  Macevoy  and  Avellaneda  considers  the  regime  (5  <?;  1  and  6ps  ~  1. 
It  gives  geometry- independent  estimates  for  L12,  extending  work  done  earlier  by  Helmholtz, 
Smoluchowski,  Fixman,  and  others. 

Ohad  Levy’s  work  with  Bergman  and  Kohn  addresses  some  issues  involving  physically 
nonlinear  behavior  of  composite  materials.  The  paper  with  Bergman  summarizes  part 
of  Levy’s  thesis  work,  concerning  optical  bistability.  The  paper  with  Kohn  explores  the 
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implications  of  duality  for  nonlinear  two-dimensional  systems.  This  work  extends  the  well- 
known  results  of  Keller  and  Dykhne  concerning  linear  composites.  It  focuses  on  three  Tna.in 
types  on  nonlinear  behavior:  (i)  the  weakly  nonlinear  regime;  (ii)  power-law  behavior;  and 
(iii)  dielectric  breakdown.  Besides  making  the  consequences  of  duality  explicit  in  each 
setting,  this  work  gives  restrictions  on  the  critical  exponents  and  scaling  functions  of  “dual 
pairs”  of  random  nonlinear  composites  near  a  percolation  threshold. 

Levy’s  work  with  Stroud  and  Palffy-Muhoray  models  the  dielectric  and  optical  response 
of  some  technologically  significant  heterogeneous  media.  The  project  with  Palfiy-Muhoray 
concerns  optical  response  of  polymer  dispersed  liquid  crystals  (PDLC’s)  -  droplets  of  liquid 
crystal  dispersed  in  a  polymer  matrix.  The  focus  here  is  on  nematic  droplets  in  an  isotropic 
polymer  matrix.  These  materials  are  being  considered  for  a  number  of  applications,  in¬ 
cluding  switchable  windows  and  active  matrix  projection  displays.  The  foundation  of  these 
applications  is  the  fact  that  the  dielectric  properties  of  a  PDLC  can  be  changed  by  an 
externally  applied  electric  field.  Levy  and  Palffy-Muhoray  have  given  the  first  satisfactory 
theoretical  model  for  this  effect.  They  use  a  Maxwell  Garnett  type  mean  field  theory  to 
represent  the  dielectric  responce  of  a  PDLC  in  terms  of  the  orientation  distribution  func¬ 
tion  of  the  droplets.  The  response  of  the  PDLC  is  a  nonlinear  function  of  the  applied  field, 
because  the  applied  field  influences  the  orientation  of  each  droplet,  changing  the  orientation 
distribution  function. 

2.7  Magnetic  materials 

•  R.V.  Kohn  and  F.  Otto,  “Small  surface  energy,  coarse-graining,  and  selection  of  mi¬ 
crostructure,”  Physica  D  107  (1997)  272-289. 

•  R.  Choksi  and  R.V.  Kohn,  “Bounds  on  the  micromagnetic  energy  of  a  uniaxial  ferro- 
magnet,”  Comm.  Pure  Appl.  Math.  51  (1998)  259-289. 

•  R.  Choksi,  R.V.  Kohn,  and  F.  Otto,  “Domain  branching  in  uniaxial  ferromagnets:  a 
scaling  law  for  the  minimum  energy,”  submitted  to  Comm.  Math.  Phys. 

The  work  of  Choksi,  Kohn,  and  Otto  addresses  some  static  and  dynamic  problems  of 
micromagnetics.  The  micromagnetic  energy  involves  a  competition  between  “bulk”  and 
“surface”  energy,  often  leading  to  complex  patterns  of  domains  and  walls.  It  is  fruitful  to 
do  a  sort  of  asymptotic  analysis,  treating  the  surface  energy  as  a  small  parameter.  The 
paper  by  Kohn  and  Otto  surveys  the  power  and  utility  of  this  viewpoint.  In  particular, 
it  provides  an  expository  summary  of  Otto’s  proposed  effective  equation  for  dynamics  of  a 
ferrofluid  in  a  Hele-Shaw  cell  under  an  applied  transverse  field. 

The  papers  with  Choksi  apply  micromagnetics  to  analyze  the  structure  of  magnetic 
domains  in  a  uniaxial  ferromagnet.  This  work  focusses  on  the  (physically  relevant)  case 
when  the  wall  energy  is  small.  Its  thesis  is  that  the  observed  branching  of  magnetic  domains 
is  a  feature  of  any  domain  pattern  with  relatively  small  energy.  This  is  demonstrated  by 
identifying  the  scaling  law  of  the  minimum  energy,  i.e.  by  proving  a  pair  of  matching  upper 
and  lower  bounds.  Furthermore  it  is  shown  that  a  domain  pattern  achieving  this  scaling 
law  must  have  average  width  of  order  where  L  is  the  length  of  the  magnet  in  the 
easy  direction.  Finally,  the  authors  argue  that  branching  is  required  by  considering  the 
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constrained  variational  problem  in  which  branching  is  prohibited  and  the  domain  structure 
is  invariant  in  the  easy  direction;  its  scaling  law  is  different.  The  earlier  paper  by  Choksi 
and  Kohn  addressed  a  2D  version  of  this  problem.  The  extension  to  3D,  achieved  in  the 
later  paper  with  Otto,  required  introduction  of  an  entirely  new  technique,  based  on  an 
interpolation  lemma  linking  L^,  L°°,  BV,  and  H~^. 

2.8  Defect  energy 

•  W.  Jin,  “Singular  Perturbation  and  the  Energy  of  Folds,”  Ph.D.  Thesis,  New  York 
University,  1997. 

•  R.  Choksi,  G.  Del  Piero,  I.  Fonseca,  and  D.  Owen,  “Structured  deformations,  fracture, 
and  hysteresis” ,  in  preparation. 

The  Ph.D.  thesis  of  Weimin  Jin,  advised  by  Kohn,  addresses  the  singularly  perturbed 
variational  problem  min /q(|VuP  —  1)^  +  e^|VVupda:dy  subject  to  u  =  0  and  du/dn  =  1 
at  dCl,  in  the  limit  as  e  0.  This  problem  arose  recently  in  work  by  Gioia  and  Ortiz  on 
blistering  of  compressed  thin  films.  It  is  the  simplest  example  of  a  singularly  perturbed 
variational  problem  whose  order  parameter  is  a  gradient  field,  and  whose  e  =  0  limit  prefers 
a  continuum  of  states.  Jin  has  obtained  new  results  concerning  the  asymptotic  energy,  and 
the  form  of  minimizers.  A  paper  based  on  this  work  is  being  prepared  for  publication. 

The  work  by  Choksi,  Del  Piero,  Fonseca,  and  Owen  addresses  the  theory  of  “structured 
deformations.”  This  developing  framework  permits  simultaneous  discussion  of  deformation 
on  two  different  length  scales,  macroscopic  and  microscopic.  The  new  paper  by  Choksi  et 
al.  clarifies  the  functional  analytic  foundation  of  the  theory  of  structured  deformations.  It 
also  applies  this  framework  to  some  specific  phenomena  from  crystal  plasticity  and  cohesive 
fracture:  yield,  hysteresis,  and  the  associated  dissipation  in  plastic  yield  are  shown  to  arise 
from  instabilities  at  the  microlevel;  also,  the  dichotomy  between  brittle  vs.  ductile  fracture 
is  clarified.  Some  detail  concerning  the  latter  point:  when  the  undeformed  length  of  a  bax 
is  larger  than  a  certain  critical  value,  the  model  predicts  that  small  extension  induces  no 
cradc  while  large  extension  induces  a  single  crack  of  a  known  size  (brittle  failure).  When 
the  undeformed  length  of  a  bar  is  smaller  than  the  critical  value,  the  model  predicts  that 
small  extension  induces  no  crack  while  large  extension  induces  a  single  crack  whose  size 
increases  gradually  from  zero  with  the  amount  of  extension  (ductile  failure). 

2.9  PDE’s  modeling  turbulent  behavior 

•  M.  Avellaneda,  “Homogenization  and  renormalization:  the  mathematics  of  multi¬ 
scale  random  media,”  in  Dynamical  Systems  and  Probabilistic  Methods  in  Partial 
Differential  Equations,  P.  Deift  et  al.  eds..  Lectures  in  Applied  Mathematics  Vol.  31, 
AMS,  1996. 

•  C.  Apelian,  R.  Holmes,  and  M.  Avellaneda,  “A  turbulent  transport  model:  streamline 
results  for  a  class  of  random  velocity  fields  in  the  plane,”  Comm.  Pure  Appl.  Math. 
50  (1997)  1053-1088. 
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•  M.  Vergassola  and  M.  Avellaneda,  “Scalar  transport  in  compressible  flows,”  Physica 
D  106  (1997)  148-166. 

The  1996  paper  by  Avellaneda  is  a  survey  of  his  work  with  Majda  and  others,  using 
methods  from  homogenization  and  probability  to  derive  rigorous  results  on  turbulent  ad- 
vection  and  anomalous  difiusion. 

The  work  of  Avellaneda  and  Vergassola  addresses  the  large-time,  large-scale  behavior  of 
turbulent  transport,  when  the  underlying  velocity  field  is  not  incompressible.  The  relevant 
PDE  is  Tt  +  V  •  (uT)  —  DAT  =  0,  where  u  =  u(x,  t)  is  a  specified  random  velocity  field 
with  mean  value  0.  When  u  is  incompressible  the  effective  equation  is  diffusive.  In  general, 
however,  transport  can  be  much  faster  due  to  the  presence  of  a  linear  term  in  the  effective 
equation  (“ballistic  transport”). 

The  work  of  Apelian,  Holmes,  and  Avellaneda  reports  analytical  and  numerical  results 
concerning  the  streamlines  of  a  specific  class  of  random  velocity  fields.  Conclusions  are 
drawn  concerning  the  statistical  topography  of  the  flow,  and  concerning  the  turbulent  dif¬ 
fusion  associated  with  it. 

2.10  Thin  film  growth 

•  T.  Schulze  and  R.  Kohn,  “A  geometric  model  for  coarsening  during  spiral-mode  growth 
of  thin  films,”  Physica  D,  submitted. 

Kohn’s  work  with  Choksi  and  Otto  addressed  the  structure  of  domains  in  a  setting 
where  energy  minimization  is  paramount.  The  work  of  Schulze  and  Kohn  addresses  the 
structure  of  domains  in  a  very  different  setting  -  spiral  growth  of  thin  films  -  where  the 
behavior  is  dominated  by  kinetic  effects.  A  convenient  model  system  is  the  growth  of  c-axis 
YBCO  high-temperature-superconductor  films  on  MgO  substrates.  Experimentalists  have 
observed  that  the  grain  size  coarsens  as  a  function  of  distance  from  the  substrate.  Schulze 
and  Kohn  conjecture  that  this  coarsening  is  the  result  of  a  geometrical  competition  between 
the  grains,  whereby  the  faster-growing  spirals  gradually  dominate.  Their  paper  formulates 
this  idea,  and  demonstrates  its  consequences  analytically  and  numerically,  for  the  simplest 
possible  model  in  which  corrections  due  to  surface  diffusion  and  curvature  are  ignored. 


3  Students 

We  directed  the  work  of  three  PhD  students  in  areas  associated  with  this  project. 

Reade  Ryan  completed  his  Ph.D.  in  July,  1996  under  the  direction  of  Marco  Avellaneda. 
His  thesis  was  entitled  “Large  Deviation  Analysis  of  Gaussian  Fields  and  the  Statistics 
of  Burgers’  Turbulence.”  He  is  currently  an  Assistant  Professor  in  the  Mathematics 
Department  at  Columbia  University. 

Weimin  Jin  completed  in  Ph.D.  in  August,  1997  under  the  direction  of  Robert  Kohn.  His 
thesis  was  entitled  “Singular  Perturbation  and  the  Energy  of  Folds”.  He  is  currently 
an  Assistant  Professor  in  the  Mathematics  Department  at  Indiana  University. 
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Matthew  Killough  will  complete  his  Ph.D.  thesis  in  August,  1998  under  the  direction  of 
Robert  Kohn.  His  thesis  is  entitled  “A  Diffuse  Interface  Approach  to  the  Development 
of  Microstructure  in  Martensite.”  He  will  begin  an  Assistant  Professorship  in  the 
Mathematics  Department  at  University  of  Minnesota  in  September,  1998. 

4  Publications 

In  refereed  journals: 

1.  M.  Avellaneda  and  P.  Swart,  “Caculating  the  performance  of  1-3  piezoelectric  com¬ 
posites  for  hydrophone  applications:  an  effective  medium  approach,”  J  Acoust.  Soc. 
Amer.  103  (1998)  1449-1467. 

2.  C.  Apelian,  R.  Holmes,  and  M.  Avellaneda,  “A  turbulent  transport  model:  streamline 
results  for  a  class  of  random  velocity  fields  in  the  plane,”  Comm.  Pure  Appl.  Math. 
50  (1997)  1053-1088. 

3.  K.  Bhattacharya,  R.D.  James,  and  P.J.  Swart,  “Relaxation  in  shape-memory  alloys. 
Part  I:  Mechanical  model,”  Acta  Materialia  45  (1997)  4547-4560. 

4.  K.  Bhattacharya,  R.D.  James,  and  P.J.  Swart,  “Relaxation  in  shape-memory  alloys. 
Part  II:  Thermomechanical  model  and  proposed  experiments,”  Acta  Materialia  45 
(1997)  4561-4568. 

5.  K.  Bhattacharya  and  R.V.  Kohn,  “Symmetry,  texture,  and  the  recoverable  strain  of 
shape-memory  polycrystals,”  Acta  Materialia  44  (1996)  529-542 

6.  K.  Bhattacharya  and  R.V.  Kohn,  “Elastic  energy  minimization  and  the  recoverable 
strains  of  polycrystalline  shape-memory  materials,”  Arch.  Rational  Mech.  Anal.  139 
(1997)  99-180. 

7.  K.  Bhattacharya,  R.V.  Kohn,  and  S.  Kozlov,  “Some  examples  of  nonlinear  homoge¬ 
nization  involving  nearly  degenerate  energies,”  to  appear  in  Proc.  Roy.  Soc.  London 
A 

8.  R.  Choksi  and  R.V.  Kohn,  “Bounds  on  the  micromagnetic  energy  of  a  uniaxial  ferro- 
magnet,”  Comm.  Pure  Appl.  Math.  51  (1998)  259-289. 

9.  R.  Choksi,  R.V.  Kohn,  and  F.  Otto,  “Domain  branching  in  uniaxial  ferromagnets:  a 
scaling  law  for  the  minimum  energy,”  submitted  to  Comm.  Math.  Phys. 

10.  R.V.  Kohn  and  F.  Otto,  “Small  surface  energy,  coarse-graining,  and  selection  of  mi¬ 
crostructure,”  Physica  D  107  (1997)  272-289. 

11.  R.V.  Kohn  and  O.  Levy,  “Duality  for  non-Ohmic  conductors  and  its  application  to 
composite  materials,”  J.  Stai.  Phys.  90  (1998)  159-189. 

12.  R.V.  Kohn  and  T.  Little,  “Some  model  problems  of  polycrystal  plasticity  with  deficient 
basic  crystals,”  SIAM  J.  Appl.  Math.,  in  press. 
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13.  0.  Levy  and  D.  Stroud,  “Macroscopic  disorder  and  the  metal-insulator  transition  in 
conducting  polymers,”  J.  Phys.:  Condensed  Matter  9  (1997)  L599-L605. 

14.  O.  Levy  and  D.  Stroud,  “Maxwell  Garnett  theory  for  mixtures  of  anisotropic  inclu¬ 
sions:  applications  to  conducting  polymers,”  Phys.  Rev.  B  56  (1997)  8035-8046. 

15.  W.  Macevoy  Jr.  and  M.  Avellaneda,  “Electoosmotic  coupling:  incorporating  larger 
surface  effects  with  a  new  length  scale,”  J.  Colloid  and  Interface  Science  188  (1997) 
139-149. 

16.  T.  Schulze  and  R.  Kohn,  “A  geometric  model  for  coarsening  during  spiral-mode  growth 
of  thin  films,”  Physica  D,  submitted. 

17.  M.  Vergassola  and  M.  Avellaneda,  “Scalar  transport  in  compressible  flows,”  Physica 
D  106  (1997)  148-166. 

Book  chapters  and  conference  proceedings: 

1.  M.  Avellaneda,  “Homogenization  and  renormalization:  the  mathematics  of  multi¬ 
scale  random  media,”  in  Dynamical  Systems  and  Probabilistic  Methods  in  Partial 
Differential  Equations.,  P.  Deift  et  al.  eds.,  Lectures  in  Applied  Mathematics  Vol.  31, 
AMS,  1996. 

2.  D.J.  Bergman  and  0.  Levy,  “Nonlinear  behavior  and  optical  bistability  in  compos¬ 
ite  media,”  in  Nonlinear  Optical  Materials,  IMA  Volumes  in  Mathematics  and  its 
Applications  101,  Springer- Verlag,  1998. 

3.  K.  Bhattachaxya  and  R.V.  Kohn,  “Recoverable  strains  in  shape-memory  polycrystals,” 
in  Proceedings  ICOMAT-95,  Journal  de  Physique  IV  -  Colloque  C8  (1995)  261-266. 

4.  K.  Bhattacharya,  R.V.  Kohn,  and  Y.-C.  Shu,  “The  Taylor  estimate  of  recoverable 
strains  in  shape-memory  polycrystals,”  in  proceedings  of  lUTAM  Symposium  on 
Transformation  Problems  in  Composites  and  Active  Materials,  Y.A.  Bahei-El-Din  and 
G.J.  Dvorak,  eds.,  Kluweri,  in  press. 

5.  P.M.  Girao  and  R.V.  Kohn,  “The  crystalline  algorithm  for  computing  motion  by 
curvatime,”  in  Variational  Methods  for  Discontinuous  Structures,  R.  Serapioni  and  F. 
Tomarelli  eds.,  Birkhauser,  1996,  pp.  7-18. 

6.  R.D.  James,  R.V.  Kohn,  and  T.W.  Shield,  “Modeling  of  branched  needle  microstruc¬ 
tures  at  the  edge  of  a  martensite  laminate,”  in  Proceedings  ICOMAT-95,  Journal  de 
Physique  IV  -  Colloque  C8  (1995)  253-259. 

7.  O.  Levy,  “Weakly  nonlinear  conductivity  and  flicker  noise  near  percolation”  in  Math¬ 
ematics  of  Multiscale  Materials,  IMA  Volumes  in  Mathematics  and  its  Applications 
99,  Springer- Verlag,  1998. 

Articles  in  advanced  stages  of  preparation: 
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1.  M.  Avellaneda,  L.  Berlyand,  and  J.-F.  Clouet,  “Frequency  dependent  acoustics  of 
composites  with  interfaces” 

2.  R.  Choksi,  G.  Del  Piero,  1.  Fonseca,  and  D.  Owen,  “Structured  deformations,  fracture, 
and  hysteresis” 

3.  R.V.  Kohn  and  B.  Niethammer,  “Geometrically  nonlinear  shape  memory  polycrystals 
made  from  a  two- variant  material” 

4.  O.  Levy  and  P.  Palffy-Muhoray,  “Dielectric  properties  of  polymer  dispersed  liquid 
crystals” 

5.  O.  Levy  and  P.  Palffy-Muhoray,  “Dielectric  response  of  polymer  dispersed  liquid  crys¬ 
tals  near  the  nematic-isotropic  transition” 


